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1.0  EXECUTIVE  SUMMARY 


1.1  Abstract 

An  overview  of  Bremsstrahlung  radiation,  which  is  prevalent  in  high-temperature  fusion 
plasmas,  is  given;  as  well  as  estimates  for  the  power  and  energy  flux  to  a  reflector  wall  from  a 
dense  plasma  focus  device  burning  p-nB.  From  these  calculations,  the  wall  temperature  and 
ablation  rate  can  be  calculated.  The  physics  derived  from  Inertial  Confinement  Fusion  (ICF) 
Hohlraum  use  is  applicable  in  the  case  where  the  wall  temperature  exceeds  105  K,  and  the 
reemission  flux  is  estimated  for  various  flux  levels.  Provided  a  method  could  be  found  for 
decreasing  the  incoming  flux  to  the  wall  (low  Z  material  shields,  cooling  methods),  multilayer 
structures  would  provide  the  best  method  for  reflecting  soft  and  hard  x-rays,  and  experimental 
results  from  high-energy  photon  experiments  are  provided. 

1.2  Summary 

An  investigation  of  Bremsstrahlung  radiation  has  been  performed  and  it  was  seen  that 
single  film  Hohlraum-like  cavities  are  best  suited  due  to  the  high  energy  flux  levels.  If  the 
reflection  cavity  were  made  small  enough  for  the  flux  levels  to  reach  10  "  W/cm“,  the  x-rays 
would  be  reemitted  10-12  times,  according  to  numerical  hydrodynamic  models  designed  at 
Sandia  National  Laboratories.  Multilayer  structures  provide  a  much  greater  range  of  energies 
over  which  photons  can  be  deflected;  however,  they  could  not  withstand  the  flux  levels  involved 
(~  10  W/cm  ),  have  not  been  designed  for  the  photon  energies  of  interest  (>  200  keV),  and  only 
reflect  at  very  small  angles  of  incidence  (~  mrad).  If  these  physics  could  be  resolved,  multilayer 
structures  would  provide  an  excellent  option  for  the  reflection  of  x-rays.  Current  tests  are  at  the 
100  keV  energy  range,  with  reflectivites  close  to  30%,  which  is  lower  than  the  50%  assumed  in 
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the  prior  dense  plasma  focus  (DPF)  study.  Further  investigation  of  inverse-Bremsstrahlung  is 
necessary  to  see  if  the  10  reemissions  found  in  Hohlraums  corresponds  to  at  least  50%  reflection. 

1.3  Introduction 

Bremsstrahlung  radiation,  which  is  German  for  “braking  radiation,”  occurs  when  charged 
particles  are  decelerated  by  collisions  with  other  charged  particles  and  emit  photons.  This  form 
of  radiation  is  pervasive  in  high  temperature  plasmas  of  fusion  interest  and  constitutes  an  energy 
loss  and  cooling  mechanism  for  the  plasma.  The  problem  of  Bremsstrahlung  emission  is 
amplified  in  our  case  because  of  the  high  temperatures  achieved  by  the  p-nB  DPF  pinch 
(~  MeV),  and  the  radiation’s  Z-  dependence. 

The  purpose  of  this  report  is  to  investigate  the  reflection  physics  for  the  high-energy 
Bremsstrahlung  radiation  emission  characterized  in  p-nB  DPF  fusion.  A  literature  search  has 
been  performed  on  Hohlraum  cavities,  multilayer  reflectors,  and  several  cooling  techniques  and 
will  be  presented.  Fortunately,  most  of  the  numerical  studies  done  at  Sandia  National  Laboratory 
on  Hohlraums  are  directly  applicable  to  our  case.  The  survey  of  the  multilayer  research  shows 
that  photon  energy  levels  of  DPF  interest  have  not  yet  been  investigated;  however,  the  general 
trend  is  towards  high  energy  reflection,  and  basic  analytical  models  are  produced.  The  vast 
majority  of  chamber  cooling  work  that  has  been  done  concentrates  on  the  transport  mechanism 
of  convection  and  ignores  radiation.  This  obviously  is  not  valid  in  our  case;  however,  the  basic 
analytical  relations  are  offered  to  give  insight  into  the  areas  needing  to  be  investigated  for 
successful  cooling.  Also,  the  mass  flow  rate  is  estimated  which  will  give  an  idea  of  how  much 
extra  propellant  needs  to  be  carried.  Finally,  the  results  found  here  are  compared  to  those 
assumed  in  the  earlier  DPF  investigation  burning  p-1  *B  [1]. 
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1.4  Bremsstrahlung  Reflection 

The  issue  of  Bremsstrahlung  reflection  is  critical  for  successful  operation  of  the  DPF 
device  [1].  The  underlying  physics  depend  on  the  energy  and  flux  of  the  incoming  radiation  to 
the  wall.  If  the  incoming  angles  and  fluxes  are  small,  then  multilayer  structures  would  provide 
the  best  option  for  reflection,  since  both  hard  and  soft  x-rays  have  been  detected  in  experiments. 
If  the  fluxes  and  wall  temperatures  reach  high  enough  levels  (>105  K),  then  a  plasma  will  be 
created  at  a  wall,  and  numerical  radiative  hydrodynamic  analysis  is  necessary  to  find  the 
reemission  of  the  photons.  Both  cases  will  be  presented  here,  beginning  with  an  estimation  of 
the  flux,  temperature,  and  ablation  at  the  wall  using  classical  heat  transfer  relations. 


1.5  Bremsstrahlung  Radiation 

In  order  to  correctly  ascertain  the  relevant  reflection  physics,  it  is  first  necessary  to  define 
the  incoming  energy  and  power  flux  to  the  wall.  The  Bremsstrahlung  emission  spectrum  is  given 
by  [2]: 


dF l 

dA 


=  6.01x10  36 


gneZeff  f- 12.40  "1 
A24kTe  {  AkTe  J 


(1) 


where  A  is  the  photon  wavelength,  ne  is  the  electron  temperature,  Zejy  is  the  effective  atomic 
number,  and  kTe  is  the  electron  temperature.  The  “Gaunt  factor”  g,  which  takes  into  account 
quantum  effects,  approaches  2n  (3)  ~  at  high  plasma  temperature  (>  550,000  K)  [3];  and  this 
value  is  used  in  the  analysis.  The  peak  Bremsstrahlung  emission  wavelength  (angstroms)  is 


6200 
kTe  (eV) 


(2) 


The  total  Bremsstrahlung  power  density  is  estimated  using  [2]: 
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(3) 


P  =5.35x10  37 Z  ~nhkT  )2, 
tot  ejje  eJ 

where  Zeff  is  defined  as: 

z»;=-Z».z.!  <4> 

ne 

and  has  a  value  of  13  for  the  p-1  *B  reaction.  The  electron  density  is  on  the  order  of  1025  m'3,  and 
the  ion  temperature  1.5  MeV,  which  gives  a  Bremsstrahlung  power  density  of  5  x  10  W/cm  . 
The  wall  loading  (input  flux)  is  set  equal  to  the  Bremsstrahlung  output  power  divided  by  the  area 
of  the  wall  that  is  exposed  to  the  plasma.  The  radiation  energy  after  each  pulse  is  found  by: 


“  TPPtotV pin 


(5) 


where  vp  is  the  pinch  lifetime,  and  Vpin  is  the  volume  of  the  pinch  region  of  the  DPF  device.  The 
power  deposited  to  the  wall  is  found  by  multiplying  the  energy  with  the  repetition  rate  v: 


=  Eu  v 
or 


(6) 


Using  the  values  found  from  the  DPF  study  [1]  for  a  500  kN  propulsion  unit,  an  energy  and 

5  2  6  2 

power  flux  of  1.23  x  10  J/cm  and  1.26  x  10  W/cm  are  found,  respectively.  These  values  will 
be  used  for  the  rest  of  the  report. 


1.6  Wall  Heating  and  Evaporation 

In  this  section  the  evaporation  rate  and  wall  temperature  will  be  estimated  using  classical 
heat  transfer.  Once  an  estimation  of  these  parameters  is  found,  a  more  accurate  analysis  can  be 
obtained.  The  analysis  will  closely  follow  the  work  done  by  Kammash  [4],  who  formulated  the 
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problem  under  conditions  of  thennonuclear  interest.  When  a  solid  is  heated  to  a  high 
temperature,  some  of  the  atoms  that  are  in  the  high-energy  tail  of  the  thermal  distribution  will 
have  sufficient  energy  to  overcome  the  surface  binding  energy;  if  the  momentum  of  these  atoms 
at  the  surface  is  directed  away  from  the  surface,  they  will  evaporate.  The  rate  of  evaporation 
can  be  estimated  from  the  vapor  pressure  of  the  solid  material  above  the  surface  and  is  given  by 
[4]: 


dn  ~  Nv 

—  =  a - 

dt  4 


(7) 


where  a  is  the  probability  that  an  atom  from  the  gas  phase  sticks  at  the  surface  (generally  set 
equal  to  1),  and  v  is  the  mean  velocity.  Noting  that  P  =  NkT,  and  that  the  mean  velocity  is 

1  /9  • 

proportional  to  (T/M)  ,  Equation  (7)  can  be  put  in  the  fonn 


dn  _  3.5x  1 022  P(T)  [Torr] 
dt  VM  Vt  (K) 


where  M  is  the  atomic  mass  number.  The  most  severe  thermal  loading  on  the  wall  occurs  when 
the  plasma  is  at  the  end  of  a  discharge  and  the  plasma  is  “dumped”  on  the  wall  in  a  very  short 
period  of  time.  In  order  to  calculate  the  number  of  atoms  evaporated  during  one  heat  pulse,  we 
must  first  calculate  the  temperature  increase  in  the  solid  as  a  result  of  sudden  heating.  This 
temperature  will  also  be  used  to  find  the  amount  of  cooling  necessary.  An  outline  and  solution 
of  the  problem  is  given,  and  the  reader  is  referred  to  Kammash  [4]  for  the  complete  analysis. 

The  equation  of  interest  in  this  case  is  the  unsteady  heat  conduction  equation  which  is 
applied  to  a  semi-infinite  solid  so  that  the  temperature  is  a  function  of  one  spatial  dimension  only 
(neglects  curvature);  i.e., 
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where 


a  = 


k 

~PCP 


(10) 


is  the  thermal  diffusivity,  k  the  thermal  conductivity,  Cp  is  the  specific  heat  at  constant  pressure, 
and  p  is  the  density  of  the  solid.  An  illustration  of  1-D  heat  conduction  and  the  associated 
boundary  conditions  are  shown  in  Figure  1. 


&T 

dx 


7  =  Oatf  =0 


7  =  0  at  x  =  oo 


Figure  1:  1-D  Heat  Conduction  Illustration 


The  boundary  conditions  are: 


T  =  0  at  t  =  0 


f)T  F 

—  =  -C  = - kA  atx  =  0  (11) 

dx  x 


T  =  0  at  x  =  oo 


The  second  boundary  condition  represents  the  heat  flow  into  the  solid  whereby  a  total  energy  Eb,- 
is  dumped  in  the  area  A  in  the  time  r.  The  complete  solution  of  Equation  (9)  during  the  heating 
phase  (t  <  x)  can  be  written  as 


T(x,t )  =  Cyfa 


W 


exp 


4  at 


erfc 


2  4at 


(12) 


where  the  function  erfc(x)  is  related  to  the  familiar  error  function  erf(x),  by 

erfc(x)  +  erf(x)  =  erf(cc)  =1  (13) 

At  the  surface  (x  =  0),  Equation  (12)  reduces  to 


F  4  r- 

mo=— ^ 

Ax  V  npkC p 


(14) 


The  maximum  temperature  occurs  at  t  =  r,  which  is  due  to  the  dumping  of  energy  per  unit  area  in 
time  t. 


F  14 

^(0,  x)  =  (AT)max  =  —  —j=. 


A  4x\  npkC, 


(15) 
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This  temperature  increase  is  of  interest  because  it  is  indicative  of  the  other  reflection  physics  at 
the  wall.  Also  of  special  interest  is  the  surface  temperature  shortly  after  the  end  of  heating, 
which  can  be  approximated  as 


T(0J)  =  Tmax  exp 


(16) 


To  find  the  number  of  atoms  evaporated  during  each  heat  pulse,  it  is  necessary  to  assume 
uniform  deposition  during  the  dumping  time  r,  constant  heat  of  sublimation  AH,  constant  thermal 
conductivity  k,  and  a  constant  specific  heat  Cp  in  the  temperature  range  of  interest.  The  number 
of  atoms  evaporated  during  each  pulse  is  obtained  by  substituting  Equation  (7)  along  with 
Equation  (8)  into  Equation  (12)  and  integrating  over  the  heating  period.  The  corresponding 
number  during  the  cooling  phase  is  obtained  using  a  similar  process.  The  end  result  is 


n  =  0.2r 


dn 

dt 


(Tmax  )  ’ 


(17) 


where  the  dn/dt  term  is  the  time  rate  change  in  atoms  evaluated  at  Tmax.  This  equation  is  used  to 
find  the  thickness  of  ablated  material  per  pulse,  which  is  shown  in  Figure  2.  The  number  of 
pulses  goes  up  to  roughly  8.5  x  105  pulses,  which  is  how  many  times  it  would  fire  if  it  ran  at  a 
repetition  rate  of  10  Hz  continuously  for  one  day. 
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Figure  2:  Ablated  Material  as  a  Function  of  DPF  Pulses 


Also  of  interest  is  the  wall  temperature  as  a  function  of  wall  “dumping”  time,  which  is  shown  in 
Figure  3. 
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Figure  3:  Log  Scale  of  Deposition  Time  vs.  Temperature 
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It  can  be  seen  that  the  temperatures  reach  above  105  K,  and  when  the  temperature  of  the  wall 
begins  to  exceed  this  value,  the  wall  will  itself  become  an  intense  radiator  and  eventually 
determine  the  radiation  field  in  the  cavity.  The  requisite  physics  have  been  studied  in  the  use  of 
Hohlraums,  which  is  discussed  in  the  following  chapter. 
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2.0  X-RAY  CONFINEMENT  PHYSICS 


2.1  Hohlraum  Cavities 

The  investigation  of  x-ray  confinement  in  fusion  systems  has  been  driven  by  research  on 
the  use  of  Hohlraums  in  Inertial  Confinement  Fusion  (ICF).  In  ICF,  small,  spherical  fuel  pellets 
are  imploded  by  high-power  lasers  or  ion  beams.  In  order  to  achieve  symmetric  compression, 
the  pellets  can  be  placed  in  cylindrical  gold-plated  cavities  called  Hohlraums.  The  Hohlraum 
contains  small  holes  through  which  beams  pass,  and  when  targeted  by  a  laser  or  ions,  the 
Hohlraum  converts  the  beam  into  soft  x-rays  on  the  inner  wall,  which  subsequently  provide 
indirect  heating  of  the  total  inner  wall.  The  confinement  effect  arises  because  the  cavity  wall 
heats  up  due  to  the  heating  from  this  source  and  becomes  itself  a  strong  emitter  of  thermal  soft  x- 
ray  radiation  [5].  In  this  way  a  fraction  of  the  flux  which  the  wall  receives  from  the  source  is 
reemitted  from  the  cavity.  A  diagram  of  the  ICF  concept  utilizing  a  Hohlraum  is  illustrated  in 
Figure  4  [6]. 


(a)  Indirect  drive 


Capsule  with 
ablator  and  fuel 


Hohlraum 

Figure  4:  Indirect  Drive  ICF  Concept 


Thus  far,  ignition  experiments  have  used  indirect  drives  in  which  an  external  laser  light  heats  the 
Hohlraum  cavity.  The  light  is  converted  with  close  to  100  percent  efficiency  into  an  intense  flux 
of  x-rays  of  almost  1,000  terawatts  per  square  centimeter  [6].  The  x-rays  converge  on  the 
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capsule's  outer  ablator  layer,  heating  and  expanding  it.  The  rocket-like  blow-off  of  the  ablator 
then  pushes  the  rest  of  the  capsule  inward,  compressing  the  interior  fuel  to  extreme  pressures  and 
temperatures. 

2.2  Numerical  Models 

Quantitative  modeling  of  the  injection  of  several  laser  beams  into  the  cavity,  the 
subsequent  conversion  of  the  laser  light  into  soft  x-rays,  and  the  resulting  spatial  distribution  of 
energy  deposition  possess  considerable  difficulties,  owing  to  the  complicated  geometry  and  the 
rather  involved  physics  of  laser  light  conversion  [5],  In  order  to  reach  a  more  tractable  model, 
the  laser  is  replaced  by  a  fictitious  source  of  x-rays  located  inside  the  cavity  wall.  In  our  case, 
this  x-ray  source  is  the  Bremsstrahlung  radiation  produced  from  the  DPF  device,  and  the  pinch 
region  will  serve  as  the  target  region  as  opposed  to  a  frozen  D-T  pellet. 

In  the  physics  model  it  is  also  assumed  that  radiation  and  matter  are  in  complete 
thennodynamic  equilibrium  in  the  cavity;  i.e.,  the  wall  emits  Planck  radiation  according  to 
Boltzmann’s  equation  into  the  cavity,  and  the  loss  of  energy  by  diffusion  into  the  wall  can  be 
calculated  by  an  approximation  of  radiation  heat  conduction.  The  reemission  of  the  x-rays  is 
detennined  by  a  nonlinear  heat  wave  which  forms  on  the  inside  of  the  wall.  A  diagram  of  the 
wave  propagation  process  is  illustrated  in  Figure  5  [7].  At  time  t  =  0  (a)  the  body  is  brought  into 
contact  with  a  thermal  bath.  For  t  >  0  (b)  first  a  nonlinear  heat  wave  runs  into  the  undisturbed 
material.  Subsequently,  hydrodynamic  motion  of  the  heated  material  becomes  important;  and 
the  heat  wave  is  overtaken  by  a  shock  wave  and  the  ablative  heat  wave  forms. 


12 


Figure  5:  Heating  of  a  Solid  Body  in  Contact  with  a  Thermal  Bath 

Additionally,  cases  have  been  studied  of  a  cavity  with  (open)  and  without  (closed)  holes  [6];  the 
latter  being  where  our  interest  lies,  since  we  won’t  encounter  external  laser  heating. 

2.2.1  Similarity  Analytical  Model 

There  have  been  varying  levels  of  sophistication  in  the  x-ray  reflection  models  produced, 
but  in  general,  the  following  is  true.  The  key  parameter  for  describing  radiation  confinement  in  a 
cavity  is  the  reemission  coefficient  of  the  x-ray  heated  wall  [5].  It  is  determined  by  a  radiation- 
driven  ablative  heat  wave  propagating  into  the  depth  of  the  wall  material,  as  described  by  Pakula 
and  Sigel  [7].  To  understand  the  situation,  consider  the  case  where  a  solid  gold  wall  is  irradiated 
from  the  left  with  a  constant  radiation  flux,  as  shown  in  Figure  6. 
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Figure  6:  Flux  at  the  Interface  Between  Wall  and  Inner  Cavity 


The  wall  receives  a  flux  Ss  from  a  source  and  an  incident  flux  S)  of  thennal  radiation  from  the 
other  wall  elements  in  the  cavity.  The  wall  radiates  a  reemitted  flux  Sr  into  the  cavity,  whereas  a 
net  heat  flux  S/,w  of  radiation  diffuses  into  the  wall  by  a  process  known  as  photon  diffusion.  The 
energy  balance  of  this  process  is  given  by: 


Ss  +  Si  =  Sr  +  Shw  (18) 

For  a  completely  closed  cavity  (no  holes),  the  source  flux  must  flow  into  the  wall;  there  is  no 
other  loss  than  the  heat  propagation  into  the  wall.  In  this  case,  the  temperature  at  the  boundary 
between  the  ablation  heat  wave  and  the  vacuum  is  given  by  the  self  similar  solutions  for  the 
ablative  heat  wave  in  gold  as  [5]: 


r  =  3.11xl06Ss4/V/13  (K)  (19) 

14  2 

where  the  units  of  flux  are  in  10  W/cm  and  the  units  of  time  in  nanoseconds.  A  plot  of  the 
wall  temperature  in  the  time  frame  of  interest  is  shown  in  Figure  7. 
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Figure  7:  Wall  Temperature  as  a  Function  of  Heating  Time 


For  blackbody  radiation,  the  reemitted  flux  .S,- has  the  fonn: 

Sl  =  art  (20) 

Therefore,  the  reflected  flux  and  temperature  can  be  obtained  from  the  source  flux  and  time.  The 
ratio  N  =  Sr  /  S/nv  is  called  the  reemission  factor  of  the  wall  and  is  a  measure  for  the  quality  of 
radiation  confinement.  The  factor  N  characterizes  the  wall  and  depends  on  the  material  as  well 
as  time.  Physically  it  corresponds  to  the  number  of  reemissions  of  the  source  energy  inside  the 
cavity.  A  plot  of  the  reemission  coefficient  with  the  input  Bremsstrahlung  flux  as  a  parameter  is 
shown  in  Figure  8. 
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Figure  8:  Reemission  Coefficient  vs.  Time  for  Various  Bremsstrahlung  Fluxes  on  Gold 

(Similarity  Model) 


It  should  be  noted  from  this  plot  that  the  number  of  reemissions  increases  strongly  with  absorbed 
radiation  flux  and  time.  A  more  complete  model  that  takes  into  account  the  atomic  structure  and 
material  opacity  is  described  below. 

2.2.2  Murakami  Computational  Model 

More  advanced  models  have  been  used  to  estimate  the  reemitted  x-ray  flux.  Again  using 
an  ablative  heat  wave  to  model  the  x-ray  flux  inside  the  reflector,  Murakami  [8]  assumed  the 
coefficient  for  radiation  thermal  conductivity  is  directly  related  to  the  radiation  mean  free  path 
and,  therefore,  to  the  Rosseland  mean  opacity.  Taking  the  power  law  approximation  for  the 
frequency  averaged  mean  free  path  (Rosseland  mean): 


£  = 


£0T 


a 


(21) 


With  temperature  T,  density  p,  and  parameters  1 0 ,  a’,  (1  Murakami  obtains  the  scaling  law  for 
the  reemitted  flux  in  the  form: 
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(22) 


S  =M  Svtg 
r  o  a 

where  Sa  is  the  absorbed  flux  and  is  set  equal  to  the  energy  flux  calculated  in  Equation  5,  and  t  is 
time.  The  constants  Ma,  v,  and  time  t,  can  be  analytically  derived  by  first  adjusting  (21)  to 
corresponding  opacity  properties,  and  then  detennining  the  reemitted  flux  in  (22)  by  using  the 
analytic  relations  given  in  Reference  9.  A  best  fit  of  these  scaling  laws  to  actual  numerical 
simulations  that  take  into  account  the  full  complexity  of  the  equation  of  state  and  opacity  tables 
in  the  appropriate  temperature  and  density  regions  has  been  tabulated  [7],  For  gold,  Equation 
(22)  becomes: 

5  =13.01SL05t0-46,  (23) 

r  a 

where  the  fluxes  are  in  units  of  10  W/cm“  and  time  is  in  the  units  of  10'  s.  An  important  result 
not  shown  here  from  Murakami’s  analysis  is  that  high  Z  gold  reemits  incident  radiation  ten  times 
more  efficiently  than  low  Z-carbon;  and  vice  versa,  carbon  absorbs  ten  times  more  per  unit  area 
than  gold  when  in  contact  with  the  same  radiation  field.  By  using  (23),  a  plot  of  the  reemission 
coefficient  vs.  time  can  be  reproduced. 
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Figure  9:  Reemission  Coefficient  vs.  Time  for  Various  Bremsstrahlung  Fluxes  For  Gold 

(Murakami  Model) 

The  flux  and  time  values  in  this  figure  are  the  same  as  in  Figure  8.  The  immediate 
difference  that  can  be  seen  from  Figures  8  and  9  is  that  the  reemission  coefficient  N  depends 
very  weakly  on  the  absorbed  flux  and  on  time  approximately  at  t  .  This  gives  N  values  a  factor 
of  2-3  times  smaller  than  the  Similarity  Model.  Experiments  have  shown  that  the  scaling  laws 
described  by  Murakami  are  in  better  agreement  with  x-ray  confinement  tests  perfonned  in  the 
laboratory,  although  improved  opacity  calculations  are  desired  [8].  Another  trend  that  is  evident 
is  that  the  reemission  coefficient  increases  with  incoming  flux,  which  would  suggest  a  smaller 
surface  area  reflector  would  be  desirable.  For  instance,  in  the  case  of  absorbed  flux  of  10  ~ 
W/cm  ,  the  x-rays  would  be  reemitted  roughly  ten  times  before  being  lost. 

There  are  certain  limitations  to  using  a  Hohlraum-like  cavity.  The  design  and  analysis 
assume  a  constant  input  flux  and  radiation  spectrum,  which  will  not  be  the  case  during  the 
collapse  of  the  pinch.  Furthennore,  it  can  be  seen  from  both  models  that  the  reemission 
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coefficient  increases  with  the  incoming  wall  flux,  which  would  suggest  the  use  of  small 
confinement  cavities.  Although  the  ablation  material  depth  is  typically  on  the  order  of  microns, 
for  longer  missions  (months  -  years)  this  can  add  up  over  hundreds  of  thousands  of  pulses 
spanning  over  an  extended  operating  period.  Cooling  of  the  Hohlraum  material  would  only 
decrease  its  performance,  as  it  would  decrease  the  plasma  intensity  inside  the  wall  and  limit  the 
number  of  reemissions.  Because  of  these  shortcomings,  the  use  of  multilayer  materials  needs  to 
be  investigated,  as  well  as  the  cooling  techniques  necessary  for  their  proper  use.  Additionally,  a 
simple  analytical  analysis  will  show  that  multilayer  structures  have  superior  reflectivities  than 
those  of  gold  alone. 
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3.0  MULTILAYER  REFLECTORS 


3.1  Introduction 

Broad-band  optics  for  x-rays  have  traditionally  consisted  of  high  Z,  high  density,  single 
element  thin- (il ms  (e.g.,  gold),  reflecting  in  grazing  incidence  by  total  reflection  [10].  A 
promising  alternative  is  the  “super-mirror  multilayer  structure,”  in  which  the  layer  spacing  has 
been  gradually  decreased  as  a  function  of  depth.  Lead,  Tungsten,  and  numerous  carbides  are  the 
materials  that  are  have  been  tested  for  these  structures.  The  lower  energy  z-rays  will  be  totally 
reflected  from  the  surface  layer,  while  the  harder  x-rays  will  penetrate  into  the  multilayer  until  a 
region  is  reached  where  the  layer  spacing  is  such  that  the  x-ray  is  reflected.  The  principle  is 
schematically  illustrated  below  [10]. 


Figure  10.  Principle  of  Broad  Band  X-ray  Reflection 


Depth-graded  x-ray  multilayers  consist  of  bilayers  comprising  material  pairs  selected  for  both 
their  optical  and  material  properties,  with  a  range  of  bilayer  thicknesses  chosen  so  as  to  reflect 
over  a  wide  energy  band;  multilayers  designed  for  use  above  100  keV  in  particular  contain 
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hundreds  to  thousands  of  nanometer-scale  layers  having  near  perfect  interfaces  in  order  to 
achieve  optimal  performance  [11]. 


3.2  Reflector  Theory 

Multilayer  structures  work  on  the  principle  of  total  reflection,  which  is  illustrated  in 
Figure  1 1  [12]: 


lll=l 


II- 


Figure  11:  The  Principle  of  Total  Reflection 


The  reflective  index  n  is  a  ratio  of  the  corresponding  perpendicular  wave  vector  components  kZ2 
and  kzi  when  a  photon  passes  an  interface  between  vacuum  and  matter.  Consider  a  photon 
traveling  within  a  vacuum  (ni  =  1)  with  the  sample  material  having  a  refractive  index  m.  In  case 
ni  is  larger  than  n2,  the  photons  will  be  refracted  away  from  the  normal  direction  (black  arrow). 
Inside  the  material,  the  perpendicular  component  of  the  photon  vector  is  reduced  to  kz2  compared 
to  its  value  kzi  outside  the  material.  As  the  angle  of  the  incident  photon  becomes  smaller,  a  so- 
called  evanescent  wave  traveling  along  the  surface  is  created  (red  arrow).  A  further  decrease  of 
the  angle  of  incidence  leads  to  total  reflection.  The  incoming  radiation  practically  cannot 
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penetrate  the  material  (green  arrow).  The  angle  of  incidence,  at  which  total  reflection  occurs,  is 


called  the  critical  angle  9C. 

Total  reflection  is  limited  to  a  regime  in  which  the  reflected  photon  momentum  transfers 
q  is  less  than  the  critical  momentum  transfer  (below  which  total  reflection  occurs)  [10]: 

pZ' 


qc  =  0.0292 


V  A  j 


(24) 


3 

where  Z  is  the  atomic  number,  A  is  the  atomic  weight,  p  is  the  density  (g/cm  )  and  units  of  qc  are 
in  A'1.  The  corresponding  critical  angles  and  energies  can  be  found  by  noting  that  for  specularly 
reflected  x-rays,  the  momentum  transfer  is 


q  = 


An  sin  0 

I 


(25) 


And  the  energy  of  incoming  x-ray  (keV)  is  given  as  a  function  of  the  wavelength  (angstroms)  by: 


12.398  ^  An 
X  ~  X 


(26) 


By  combining  (25)  and  (26)  one  obtains  q  =  E  sin  6.  This  implies  that  high  photon  energy 
applications  require  reflection  at  very  small  angles,  which  in  turn  requires  very  long  mirrors  with 
very  stringent  figure  criteria.  A  gold  mirror  reflecting  up  to  80  keV  would  thus  require  a  1  mrad 
angle  of  incidence  and  a  length  of  1  m  to  reflect  a  1  mm  wide  beam. 

Due  to  interface  effect  and  the  proximity  of  different  materials,  multilayer  structures  have 
unique  properties  different  from  single  film  materials  [13].  Depth  graded  x-ray  multilayers 
consist  of  bilayers  comprising  material  pairs  (i.e.,  W/SiC)  selected  for  both  their  optical  and 
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material  properties,  with  a  range  of  bilayer  thicknesses  chosen  so  as  to  reflect  over  a  wide  energy 
band.  The  peak  reflectance  attainable  from  an  x-ray  multilayer  depends  in  practice  on  the 
reflection  coefficient  (ratio  of  reflected/scattered  power  to  incident  power)  at  each  interface, 
determined  by  the  optical  constants  (e.g.,  opacity)  of  the  materials  and  by  the  interface  width  that 
characterizes  the  degree  of  interface  perfection  [12].  The  energy  of  the  hardest  x-ray  that  may  be 
reflected  by  the  structure  can  be  found  by  utilizing  Bragg’s  Law: 

^  =  Xnin  sin  0  (27) 

in  which  dmi„  is  the  thickness  of  the  bottom  layer.  By  combining  this  with  the  relation  given  by 
(26)  we  are  able  to  obtain 


^-?-  =  Esm0xq  (28) 

d 

min 

This  value  of  q  is  almost  four  times  higher  than  the  value  of  a  gold  film  reflector.  The  reflection 
band  of  an  optimal  supermirror  may  therefore  be  expected  to  be  a  factor  of  4  wider  than  the  band 
of  a  thin-film  reflector.  However,  this  additional  bandwidth  is  gained  at  the  cost  of  reflectivity; 
absorption  reduces  the  reflectivity  at  energies  where  reflection  happens  within  the  supermirror. 
The  added  bandwidth  outweighs  the  lost  reflectivity  in  the  case  of  Bremsstrahlung  radiation, 
since  both  soft  and  hard  x-rays  are  present. 

3.3  Experimental  Findings 

To  illustrate  the  advantages  of  multilayers  over  single  film  reflectors,  a  comparison  has 
been  done  measuring  the  20  -  95  reflectivity  of  an  Au-coated  reflector  and  a  600  element  bi¬ 
layer  W/Si  supennirror.  The  results  of  the  study  are  shown  in  Figure  12  [10]. 
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Figure  12:  Comparison  of  Reflectivities  of  W/Si  Supermirror  and  Gold  Film 


In  comparing  the  measured  reflectivities  of  the  supermirror  and  the  gold  coating,  the 
ability  of  the  supermirror  to  reflect  at  higher  energies  is  clearly  seen.  For  angles  of  4.5  and  3.0 
mrad,  the  cutoff  of  the  W/Si  multilayer  is  determined  by  the  absorption  at  69.5  keV.  At  1.5 
mrad,  the  absorption  edge  is  seen  to  have  little  effect,  and  one  may  therefore  expect  considerable 
reflectivity  (~  200  keV)  at  this  angle.  The  major  limitation  of  the  supermirror,  i.e.,  a  reflectivity 
far  from  unity,  is  also  clearly  illustrated  by  the  measurements.  This  is  especially  severe  for  the 
higher  angles,  but  at  2  mrad  and  below,  a  reflectivity  above  30%  is  obtained  in  the  whole  band. 
In  order  to  overcome  this  problem,  other  material  combinations  could  be  used  in  the  mirror,  and 
are  being  investigated.  Using  Equation  (2),  it  is  seen  that  the  highest  energy  photons  will  be 
coming  in  with  an  energy  of  over  200  keV,  which  is  above  the  limits  currently  known.  However, 
W/  SiC  supermirrors  have  already  been  produced  [12]  and  have  been  shown  to  reflect  well  above 
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100  keV.  Also,  if  the  bulk  of  x-rays  are  going  to  be  in  a  narrow  bandwith,  a  multilayer  could  be 
fabricated  with  a  larger  reflectivity  in  this  band  of  interest. 

Although  the  added  bandwith  would  be  a  major  advantage  over  single  film  reflectors, 
there  are  several  shortcomings.  The  first  problem  is  that  the  incoming  flux  from  the  DPF  pinch 
would  quickly  destroy  the  multilayer.  A  combination  of  cooling  methods  (e.g.,  regenerative  and 
film)  might  aid  in  solving  this  problem  and  are  discussed  in  the  next  chapter.  However,  it  is 
likely  that  additional  approaches  may  be  necessary,  such  as  shield  absorber  made  of  ceramics  or 
other  low  Z  material.  Also  the  photon  energies  of  interest  (>  200  keV)  have  not  yet  been 
examined.  The  final  major  concern  is  the  small  angles  (~  mrad)  necessary  for  reflection.  This 
consideration  would  make  a  reflector  much  larger  than  the  pinch  region  unnecessary,  since  any 
radiation  coming  in  at  a  large  angle  would  not  be  reflected.  If  these  obstacles  can  be  overcome, 
then  the  multilayer  would  provide  an  attractive  alternative  to  single  film  layers.  The  reflectivity 
of  30%  obtained  in  current  experiments  is  lower  than  the  50%  assumed  in  the  original  DPF 
report  [1];  however,  the  advances  cited  may  increase  this  number  in  future  experiments. 
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4.0  WALL  COOLING  METHODS 


4.1  Introduction 

The  current  literature  on  multilayer  structures  deals  with  flux  rates  much  lower  that  those 
expected  with  typical  DPF  operation.  At  the  flux  level  and  wall  temperatures  involved,  the 
multilayer  would  be  quickly  destroyed;  therefore,  an  active  cooling  system  will  be  necessary. 
Much  of  the  work  that  has  been  done  in  this  area  relates  to  the  cooling  requirements  of 
combustion  chambers  in  chemical  rockets.  Concepts  developed  to  cope  with  this  problem,  either 
singly  or  in  combination,  include  regenerative  cooling,  radiation  cooling,  film  or  transpiration 
cooling,  ablation,  and  inert  or  endothennic  heat  sinks  [14].  For  the  temperatures  encountered 
with  the  DPF  device,  it  will  be  necessary  to  use  a  combination  of  regenerative  cooling  and  liquid 
film  cooling.  Many  liquid  rocket  engines  employ  a  film  of  liquid  fuel  as  thermal  protection  for 
the  combustion  walls. 

4.2  Film  Cooling 

The  liquid  film  cooling  process  was  experimentally  studied  in  the  1950’s  and  1960’s; 
however,  no  general  theoretical  model  was  ever  developed.  The  classical  approach  of  solving 
liquid  film  cooling  problems  is  a  turbulent  flat  plate  correlation  using  either  Reynolds’s  Analogy, 
which  assumes  identical  velocity  and  temperature  profiles  in  the  laminar  vapor  sublayer,  or 
Colburn’s  equation  based  on  the  l/7th  power  law  to  characterize  the  velocity  profile  [15]. 
However,  the  fonnulation  of  this  problem  is  difficult  to  obtain  because  of  the  complex 
phenomena  that  characterize  the  flow  of  a  high  velocity  gas  over  a  liquid  film.  Additionally, 
most  of  the  investigations  have  been  characterized  by  ambient  pressure  and  temperature, 
essentially  zero  heat  transfer  rates,  and  relatively  low  interfacial  shear  forces.  Such  conditions 
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are  significantly  different  from  those  that  characterize  the  typical  application  of  liquid-film 
cooling,  especially  at  the  DPF  conditions  considered  here.  Also,  radiative  transfer,  which  is 
prevalent  in  our  case,  is  not  for  most  chemical  rocket  chambers;  and  therefore  has  received  very 
little  attention.  Due  to  the  lack  of  experimental  and  theoretical  data  in  our  operating  regime,  only 
the  basic  relations  will  be  presented  to  give  an  order  of  magnitude  estimate  for  the  cooling 
requirements. 

In  liquid  cooled  systems,  fundamental  problems  include  the  coolant  mass  flow  rate 
required  to  cool  the  desired  internal  surface  area,  the  intact  liquid  film  length,  and  the  effect  of 
the  cooling  process  on  the  perfonnance  of  the  engine  [16].  The  basic  analytical  and  design 
problem  can  be  stated  as  that  of  determining,  for  given  liquid  and  gas  flow  parameters,  1)  the 
rate  of  coolant  injection  required  to  establish  a  desired  wetted  area;  and  2)  the  degree  of  the 
insulating  effect  of  the  gas-vapor  layer  downstream  of  the  liquid  film.  The  primary  objectives  of 
a  liquid  film  cooling  model  are  to  predict  the  temperature  profile  along  the  chamber  wall  and  to 
determine  the  film-cooled  length  to  ensure  that  sufficient  liquid  film  is  injected  across  the 
reflector.  By  looking  at  the  mass  transfer  properties,  the  necessary  mass  flow  rate  can  be 
estimated. 

4.2.1  Mass  Transfer 

The  mechanisms  involved  in  liquid  film  cooling  are  depicted  in  Figure  13  [15].  Heat  is 
transferred  from  the  hot  free-stream  gas  to  the  liquid  film  by  both  radiation  and  convection.  Heat 
energy  from  the  hot  gas  stream  increases  the  sensible  enthalpy  of  the  liquid  by  radiation, 
convection  and  conduction.  After  the  saturation  temperature  of  the  liquid  film  is  reached,  the 
incident  heat  is  used  to  vaporize  the  coolant.  The  liquid  film  tenninates  at  some  point 
downstream  of  the  injector  as  a  result  of  evaporation  and  its  entrainment  into  the  core  gas  stream. 
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The  distance  from  the  injector  to  the  tennination  point  is  the  liquid  film-cooled  length. 
Downstream  of  the  liquid  tennination  point,  the  vapor  provides  thermal  blockage  through 
gaseous  film  cooling. 


cl  condj  Cp,  T, 


Figure  13.  Control  Volume  for  Interfacial  Energy  Balance 

This  energy  is  absorbed  by  vaporizing  the  liquid  in  the  protective  film  on  the  wall.  The  vapor¬ 
generated  flow  is  known  as  the  transpiration  cooling  process.  Downstream  of  the  liquid  film,  the 
vapor  mixes  with  the  free  stream  gas  entrained  in  the  boundary  layer,  lowering  the  wall 
temperature  through  the  well-known  “gaseous  film  cooling”  process.  This  provides  thermal 
protection  downstream  of  the  dryout  point. 

The  total  heat  flux  due  to  both  convection  and  radiation,  Q,oh  is  absorbed  in  the  liquid 
film,  causing  an  initial  temperature  rise  [17]: 
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(29) 


dT«,  Q,„ 

dx  r  c  p 

where  /"is  the  local  liquid  mass  flow  rate  per  circumference.  After  the  liquid  film  reaches  the 
saturation  temperature  Tv,  the  evaporation  rate  is: 
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where  A  is  the  latent  heat  of  evaporation.  In  most  rocket  combustion  chambers,  the  radiant  heat 
flux  is  negligible  in  comparison  to  the  convective  heat  flux  [17].  However,  this  will  be  the 
complete  opposite  for  typical  DPF  operation,  and  the  radiant  heat  flux  will  largely  determine  the 
evaporation  rate  as  shown  in  Equation  (2).  The  radiant  heat  flux  can  be  estimated  through  the 
use  of  Equation  (6)  divided  by  the  surface  area  of  the  exposed  wall. 

Q  rad  =PBr!  Area-  107W/cm2  (31) 

This  flux  level  is  well  above  levels  found  in  chemical  rockets.  Because  of  this  tremendously 
high  flux,  another  concern  is  that  some  of  the  radiant  heat  can  penetrate  the  liquid  film  and  be 
absorbed  directly  at  the  combustion  chamber  walls.  With  high  enough  heat  flux,  the  liquid  film 
can  “burnout,”  as  in  nonnal  pool  boiling.  Monde  and  Katto  [18]  have  studied  this  problem  for 
heat  fluxes  on  the  order  of  106  W/m2  and  have  correlated  the  burnout  heat  flux  as: 
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where  pv,  and  pi  are  the  densities  of  the  vapor  and  liquid,  cr  is  the  surface  tension,  A  is  the  latent 
heat  of  vaporization,  L  is  the  length  of  the  heated  surface,  and  U  is  the  average  velocity  of  the 


29 


liquid  film.  It  has  been  realized  experimentally  that  when  the  critical  heat  flux  is  exceeded,  the 
liquid  film  separates  from  the  heated  surface.  However,  it  is  possible  that  this  separation  from 
the  surface  in  the  experiments  may  have  been  due  to  the  constant  heat  flux  condition  imposed. 
In  DPF  operation  the  liquid  film  may  be  less  susceptible  to  burnout  than  those  in  the  heat  transfer 
experiments,  because  separation  of  the  liquid  film  from  the  surface  would  generate  a  region  of 
droplets  and  bubbles  which  would  cause  scattering,  decreasing  the  radiative  transmission  through 
the  film.  Use  of  these  burnout  correlations  is  questionable,  since  they  are  expressed  in  terms  of 
the  overall  heated  length.  Ideally,  the  burnout  point  should  be  expressed  in  terms  of  local 
conditions,  such  as  the  local  film  thickness.  This  is  why  this  analysis  is  mainly  useful  as  an  order 
of  magnitude  calculation. 

The  convective  heat  flux  can  be  estimated  using  [17]: 

Qconv  =  h0^T ,  (33) 

where  h0  is  the  local  convective  coefficient  and  can  be  calculated  from: 

h0=KtGCpgSt0,  (34) 

where  iffis  the  turbulence  correction  factor  given  by: 

Kt=  1  +  4ef  (35) 

The  parameter  e,  has  a  value  of  0.05  -  0.2  depending  on  the  turbulence  intensities  and  distances 
from  the  liquid  film  injector,  and  we  will  use  a  value  of  0. 1 .  St„  is  the  Stanton  number  defined 
as: 

St0=/2Cf  Pr-°'6  (36) 
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The  Stanton  number  is  a  dimensionless  parameter  typically  made  up  of  other,  more  familiar 
dimensionless  parameters.  It  can  be  defined  for  heat  transfer  or  for  mass  transfer.  C/is  the 
friction  factor  expressed  as: 

Cj-  =0.0592  Re"0,2  (37) 

Rex  is  the  Reynolds  number  based  upon  the  distance  x  from  the  leading  edge  and  is  valid  for 
values  >  1.10  x  107  [15]: 

Rer  =  G—  (38) 


where  G  is  the  free  stream  gas  flow  (pg  Ug  ).  In  order  to  take  into  account  growing  boundary 
layer  present  in  the  flow,  x  is  replaced  by: 


The  convective  coefficient  is  dependant  on  the  Reynolds  number  (mass  velocity).  Table  1  shows 
how  it  varies  at  a  distance  of  15  cm  (~  half  the  reflector  distance)  from  the  liquid  injector  point 
using  Equations  34  through  39  with  hydrogen  as  a  coolant.  The  Reynolds  numbers  are  provided 
up  to  1 07,  which  is  the  upper  experimental  limit  for  which  the  equations  have  been  verified. 
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Table  1.  Convective  Transfer  Coefficient  Variation  with  Reynolds  Number 


Reynolds 

Number 

Heat 

Transfer 

Coefficient 

1.00E+05 

0.05 

2.50E+05 

0.10 

5.00E+05 

0.17 

7.50E+05 

0.24 

1.00E+06 

0.30 

2.50E+06 

0.63 

5.00E+06 

1.10 

7.50E+06 

1.52 

1.00E+07 

1.91 

From  these  values  of  the  heat  transfer  coefficient,  we  can  see  that  the  lower  the  Reynolds 
number,  the  lower  the  convective  heat  flux  will  be.  If  the  incoming  radiation  is  on  the  order  of 
106  K,  then  from  Equation  33: 

Qconv  ~  105  W/cm2  (40) 


This  produces  a  liquid  evaporation  rate  (per  area)  of: 
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We  can  thus  see  the  engineering  difficulties  encountered  with  a  cooling  system  designed  for  such 
a  significant  incoming  flux;  this  is  a  significant  mass  flow  rate  for  such  a  small  area.  The  liquid 
film  thickness  and  average  velocity  can  be  estimated  using  the  laminar  “Couette  flow”  result: 


8  = 


2//T 
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(42) 


where  xm  is  the  wall  shear  stress  and  can  be  calculated  from: 


r,=  V1CfG{Us-U,) 


(43) 
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and  the  liquid  velocity  can  be  found  from: 


U 


liq 


r 

pS 


(44) 


4.2.2  Heat  Transfer 

In  order  to  estimate  the  film  cooling  length,  the  heat  transfer  relations  derived  by 
Stechman  [19]  will  be  utilized.  Stechman  predicted  the  heat  transfer  coefficient  from  the 
combustion  gas  to  the  film  coolant  and  from  the  film  coolant  to  the  wall  by  modifying  the  one¬ 
dimensional  Bartz  equation  to  take  into  account  the  effect  of  the  liquid  and  gaseous  film. 
Equation  (45)  calculates  the  heat  transfer  coefficient  from  the  film  coolant  to  the  wall: 
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where  pi  is  the  liquid  dynamic  viscosity,  D  is  the  reflector  diameter,  Pr  is  the  Prandtl  number,  A 
is  the  reflector  area,  Hr  and  Hw  are  the  recovery  and  wall  enthalpy,  7,- and  Tw  are  the  recovery  and 
wall  temperature,  and  £  is  a  parameter  which  accounts  for  the  static  property  change  as  given  by: 
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The  heat  transfer  coefficient  for  the  liquid  is  given  by: 


h,  =  0.0288 
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where  Cpi  is  the  liquid  heat  capacity,  and  x  is  the  axial  distance.  The  film  cooling  length  was 
also  estimated  by  Stechman,  as  shown  in  Equation  48.  It  can  be  seen  from  this  equation  that  the 
length  is  proportional  to  the  heat  capacity  and  the  heat  of  vaporization: 
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ilmi  Cp,(TSqt  -  Tj  )  Omi  II v 

PKVr-Tsat)  Phg(Tr  ~Tsat)  ’ 

where  Q  is  an  empirical  correction  factor  with  a  value  between  0.5  -  1,  P  is  the  chamber 
perimeter,  Tsat  is  the  saturation  temperature,  and  Hv  is  the  total  enthalpy  of  the  vapor.  The  first 
tenn  on  the  right  represents  the  distance  from  the  injection  point  for  which  the  effective  gas 
temperature  varies  from  its  initial  injection  temperature  to  its  saturated  liquid  temperature,  and 
the  second  tenn  represents  the  length  which  is  required  for  the  hot  gas  to  completely  vaporize  the 
film.  Using  the  appropriate  values  for  the  parameters  evaluated  in  Equation  (48),  a  film  cooling 
length  of  3.2  cm  is  calculated,  assuming  a  liquid  mass  flow  rate  of  15  kg/s. 

4.3  Regenerative  Cooling 

Regenerative  cooling  is  the  process  where  a  coolant  is  passed  through  a  channel  adjacent 
to  a  wall  to  lower  the  wall  temperature.  The  steady  state  heat  transfer  through  the  chamber  wall 
of  a  liquid-cooled  rocket  chamber  can  be  treated  as  a  series  type,  steady-state  heat  transfer 
problem  with  a  large  temperature  gradient  across  the  gaseous  film  on  the  inside  of  the  chamber 
wall,  a  temperature  drop  across  the  wall,  and  a  third  temperature  drop  across  the  film  of  the 
moving  cooling  fluid.  It  is  a  combination  of  convection  at  the  boundaries  of  the  flowing  fluids 
and  conduction  through  the  chamber  walls.  The  general  steady-state  heat  transfer  equations  for 
regeneratively  cooled  thrust  chambers  can  be  expressed  as  [20]: 


q  =  h(Ta  -  Tj)  =  Q/ A 


1  lhg  +tw/K  +  \/h, 

=  hg(Ta-Twg)  ,  (49) 

=  (*/tJ(Twg  -Twl) 

=  kl<T«-Tl) 
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where  q  is  the  heat  transferred  per  unit  area  per  unit  time,  Tg  the  absolute  chamber  gas 
temperature,  7}  the  absolute  coolant  liquid  temperature,  Tw/  the  absolute  wall  temperature  on  the 
liquid  side  of  the  wall,  Twg  the  absolute  wall  temperature  on  the  gas  side  of  the  wall,  h  the  overall 
film  coefficient,  tM,the  thickness  of  the  chamber  wall,  and  k  the  conductivity  of  the  wall  material. 
An  illustration  of  these  parameters  is  shown  in  a  qualitative  temperature  profile  diagram  in 
Figure  14  [20]. 


Figure  14:  Qualitative  Temperature  Profile  of  Radiation  Reflector 

The  film  coefficients  can  be  calculated  through  the  use  of  Equations  45  and  47.  The  important 
quantities  for  controlling  the  heat  transfer  across  a  rocket  chamber  wall  are  the  fluid  film 
boundaries  established  by  the  combustion  products  on  one  side  of  the  wall  and  the  coolant  flow 
on  the  other.  The  gas  film  coefficient  largely  detennines  the  numerical  value  of  the  heat  transfer 
rate,  and  the  liquid  film  largely  determines  the  value  of  the  wall  temperatures.  The  necessary 
regenerative  cooling  is  heavily  dependant  on  the  amount  of  heat  the  liquid  film  on  the  surface  is 
able  to  take  away. 
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4.4  Limitations  of  Cooling  Models 


The  major  limitation  to  the  liquid  cooling  models  developed  and  used  is  that  there  is  no 
theory  or  experimental  work  done  in  our  area  of  interest;  chemical  rocket  chamber  conditions  are 
vastly  different  from  fusion  chamber  conditions.  It  is  therefore  unclear  how  relevant  the 
equations  presented  are.  Additionally,  the  convective  transfer  equations  are  all  derived  with  the 
assumption  that  radiation  is  negligible,  which  very  well  might  prove  to  be  the  complete  opposite. 
However,  the  analysis  presented  may  provide  order  of  magnitude  estimates  for  the  cooling 
requirements,  in  which  case  it  can  be  seen  that  a  very  high  injection  speed  and  mass  flow  rate  (30 
kg  s'1  cm'2  )  will  be  necessary  to  counteract  the  extremely  high  temperatures  found  near  the  wall 
regions. 


36 


5.0  CONCLUSIONS 


For  a  500  kN,  2000  Isp  dense  plasma  focus  propulsion  device,  the  energy  and  power  are 
1.23  x  10“  J/cnr  and  1.26  x  10  W/ctrT.  The  flux  to  the  wall  depends  on  the  exposed  area  of  the 
reflector,  although  it  has  been  shown  that  it  would  be  advantageous  to  have  a  small  reflector. 
This  is  because  the  reemission  increases  with  flux  in  the  case  of  a  gold  film  reflector;  and  in  the 
case  of  multilayers,  maximum  reflectivity  occurs  at  small  angles.  Additionally,  a  small  reflector 
would  greatly  decrease  cooling  needs.  The  use  of  gold  film,  Hohlraum-like  cavities  has  been 
explored;  and  for  an  incident  flux  of  10  W/cm  ,  the  radiation  will  be  reemitted  approximately 
10  times  before  being  lost,  according  to  the  numerical  work  done  by  Murakami.  In  the  case  of 
multilayer  structures,  x-ray  energies  of  DPF  have  not  yet  been  investigated,  although  the  trend  is 
moving  in  that  direction.  Despite  the  advantages  of  multilayer  structures  over  single  element 
layers,  such  as  a  greater  energy  band  of  reflectance,  multilayer  structures  would  not  be  able  to 
handle  the  incoming  energy  flux  levels  characteristic  of  the  p-1  *B  pinch;  and  it  is  unclear  whether 
low  Z  materials  or  cooling  methods  could  aid.  This  is  due  in  part  to  the  lack  of  theory  and 
experiment  in  cooling  methods  in  DPF-like  environments.  A  50%  reflection  rate  was  assumed  in 
the  prior  DPF  study,  which  seems  out  of  reach  of  current  multilayer  capabilities,  but  within  the 
reach  of  single-film  Hohlraum  cavities.  An  investigation  of  inverse-Bremsstrahlung  is  necessary; 
if  there  are  10  passes  of  photons  before  being  lost,  it  is  very  possible  that  the  radiation  will  be 
reabsorbed  in  the  pinch  region. 
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